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§ ■ Abstract 



The retarded van-der- Waals potential, as first obtained by Casimir and 
CN ' Polder, is usually computed on the basis of nonrelativistic QED. The hamil- 

tonian describes two infinitely heavy nuclei, charge e, separated by a dis- 
', tance R and two spinless electrons, charge — e, nonrelativistically coupled 

to the quantized radiation field. Casimir and Polder use the dipole approx- 
_C ■ imation and small coupling to the Maxwell field. We employ here the full 

& . hamiltonian and determine the asymptotic strength of the leading —R~ 7 

potential, which is valid for all e. Our computation is based on a path 
integral representation and expands in 1/R, rather than in e. 
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°2 1 Introduction 

Neutral atoms and molecules interact through the long range, attractive van der 
Waals potential which has a decay as —R~ 6 for large separation R. The quantum 
origin of this force was first recognized by London pQ. Lieb and Thirring [2] 
supply a non-perturbative proof valid for very general charge configurations. If 
one goes beyond the static Coulomb interaction in quantizing the Maxwell field, 
then the action is no longer instantaneous but travels with the speed of light 
between atoms. In a now very famous paper [3] Casimir and Polder establish 
that thereby the effective interaction potential decays somewhat faster, namely 
as — R~ 7 , which is known as the retarded van der Waals potential. For two 
hydrogen atoms, the cross over between R~ 6 and R~ 7 sets in at roughly 100 Bohr 
radii. There are both direct and indirect measurements [I] which confirm the 
theoretical prediction. 

The starting point of Casimir and Polder is nonrelativistic QED for two atoms 
separated by a distance R. Within dipole approximation they expand to fourth 
order in the coupling to the Maxwell field and obtain a prefactor of —R~ 7 which 
is proportional to the square of the electric dipole moment of a single atom. 
Later on alternative routes and simplified derivations were proposed. For an 
extensive discussion we refer to the book by Milonni [5], see also the monograph 
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by Margenau and Kestner [B] and the lecture notes of Martin and Biinzli [7J. 
Feinberg and Sucher [HUH] reconsider the issue by employing a dispersion-theoretic 
approach. Their prefactor turns out to be quadratic in the electric and magnetic 
dipole moment of a single atom. Somewhat later Boyer pU] rederived the same 
prefactor using quantum zero-point energy and semiclassical expressions for the 
level shifts due to the presence of the atoms. In our note we stick to nonrelativistic 
QED, no dipole approximation and no assumption on small coupling, and expand 
in 1/R. We use the path integral formulation, in which the subtraction of the 
ground state energy at R = oo is particularly transparent. As in previous studies 
the strength of the retarded van der Waals potential is quadratic in the electric 
and magnetic dipole moments, but with modified coefficients as compared to 

pun]. 

In mathematical physics there has been a revived interest in nonrelativistic 
QED [UJ. It is conceivable that some parts of the argument can be elevated to 
a rigorous proof. In our paper we mostly ignore this line of research, but will 
provide a more detailed discussion in the conclusions. 



2 Hamiltonian and van der Waals potential 

We consider a single hydrogen atom with an infinitely heavy nucleus located at 
the origin. The nucleus has charge e, e > 0, the electron has charge — e. We 
will use units in which h — l,c = 1, and the bare mass of the electron m — 1. 
In Section [5] we will restore the proper physical units. Let x, p be position and 
momentum of the electron. Then the nonrelativistic QED hamiltonian for this 
system reads 

H = ^(p-eA{x)) 2 -e 2 V{x) + H i . (2.1) 

For H to make sense the electron is assumed to have a prescribed charge distri- 
bution if with the following properties: ip is normalized, J dxip(x) = 1, rotation 
invariant, tp(x) = <^ ra d(|x|), of rapid decrease, and its Fourier transform, <p, is 
real. Then V is the smeared Coulomb potential 

V(x) = J dk |^(ib)| 2 |fe|- 2 e- i *- x . (2.2) 

A(x) is the quantized vector potential and iff is the field energy. These are defined 
through a two-component Bose field a(k, A), k e IR 3 , A = 1, 2, with commutation 
relation 

[a(k, A), a(k', A')*] = 5 X \>5(k - k'). (2.3) 

Explicitly 

H { =^2 [ dku{k)a{k,\)*a{k,\) (2.4) 
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with dispersion relation 

u(k) = \k\ (2.5) 

and 

A(x)=Y* [ dk e(k,X)(e ikx a(k,X)+e- ikx a(k J \y) (2.6) 

, oil 3 \/2u(k) 



A=l,2 ' 



with the standard dreibein s(k, l),s(k, 2), k = k/\k\. Thus the Hilbert space for 
H is 

H = L 2 (R 3 x )®d, (2.7) 

where $ is the bosonic Fock space over L 2 (R 3 ) £g) C 2 . From the quantization of 
the classical system of charges coupled to the Maxwell field it follows that for the 
smearing of A(x) and of V the same charge distribution has to be used. We refer 
to [H] for details. As proved by Griesemer, Lieb, and Loss [T2 l [T3 l fT^j. H has a 
unique ground state, denoted here by ip, with ground state energy E, Hip = Eip. 

The asymptotic strength k of the van der Waals potential depends on the 
properties of a single hydrogen atom only through its electric and magnetic dipole 
moment, oe and au- They are defined through the energy, W, of our system for 
weak external uniform electric and magnetic fields according to 

W = -\avEl - \a M Bl. (2.8) 

To say, H is perturbed by eE ex -x and the vector potential is perturbed by ^B cx Ax. 
Then by second order perturbation theory it follows that 



e 2 . 



a E = 2(l(ip,x-(H-E)- 1 x^)) 
«M = -^Q(^xV>))e 2 . (2.9) 

As a convention, (-, ■) denotes always the inner product on the respective Hilbert 
space. 

To investigate the van der Waals potential we consider two hydrogen atoms, 
one located at and the other at r = (0,0, R),R > 0. It will be convenient to 
define the position of the second electron relative to r. Then x\, X2+r are positions 
and pi,P2 the momenta of the two electrons. The two-electron hamiltonian reads 

H R =-{pi - eA(xi)) 2 - eV(xi) + -{p 2 - eA(x 2 + r)) 2 - e 2 V{x 2 ) 

+ H { + e 2 V R (x u x 2 ) (2.10) 
with the interaction potential 

Vr(xi, x 2 ) = —V(xi - r) - V(x 2 + r) + V(r) + V(r + x 2 - x{) 

dk \0(k)\ 2 e ik - r \k\- 2 (l - e~ ik - Xl )(l - e ifc - X2 ). (2.11) 
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H R acts on the Hilbert space L 2 (IR^ i ) ® L 2 (IR^ 2 ) <g> H R has a unique ground 
state with energy E(r). 

In the Born-Oppenheimer approximation E(r) is the effective potential be- 
tween the two neutral hydrogen atoms in their ground state. Thus the issue at 
hand is to investigate E(r) for large R. By rotation invariance E(r) depends only 
on |r| = R and we also write E(r) = E(R). For R — > oo the two atoms become 
independent and one can show that E(R) converges to 2E. The Casimir-Polder 
result is that, for small e, 

i ta o ^( BM -2 B ) = -|(^) 2 (Ia E , at ) 2 . (2.12) 

The factor (l/27r) 2 results from our use of the Lorentz-Heaviside units. We remark 
that Casimir and Polder omit in their definition of «e the factor 2, which accounts 
for the extra 1/2. «E,at is the dipole moment of a decoupled hydrogen atom. It 
is defined through 

« E ,at = 2(|(^ at ,x • (H at - E^x^e 2 (2.13) 

with 

H at = ^p 2 -e 2 V(x) (2.14) 

and ip at the ground state of the hydrogen atom, H^at = -E'atV'at- 

In our set-up, the natural dimensionless coupling constant is the Sommerfeld 
fine-structure constant 

e 2 

a = ——. (2.15) 

The energy unit is set by the ionization energy of the hydrogen atom, which is 
a 2 mc 2 and the length unit is the Bohr radius r-s = h/amc. Anticipating a decay 
as R~ 7 , the dimensionless coupling strength, k, is defined through 

E(R) -2E = -K(a, \ c A)amc 2 (R/r B y 7 (2.16) 

valid for large R. k depends on a and on the ultraviolet cutoff A in units of 
the inverse Compton wave length A c = h/mc. For this interpretation the form 
factor (p is chosen as <p(k) = </?i(A -1 &;), where (p\ rapidly interpolates between 
(p x = (2tt)- 3 / 2 for |jfe| < 1 - 6 and 1 = for \k\ > 1. 

The goal of our note is to obtain an exact expression for the strength k. 
Readers not so much interested in the details of the computation may skip ahead 
to Section [5] where the result is discussed. 

3 Path integration 

As noted by Feynman [15], in the functional integral representation of e~ tHa ,t > 
0, the interaction with the radiation field is linear in A. Therefore one can carry 
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out the Gaussian integration over the fluctuating photon field. This form will 
be particularly convenient for the Born-Oppenheimer energy E(R) — 2E. After 
such a detour we will return to operators. Our notation is formal, but rigorous 
versions are available [T6l [T7] . 

We denote by qj(t) £ M 3 the path of electron j. In case of a single electron 
we omit the index j. For the ground state energy of the two-electron system one 
obtains 



E(r) = - lirn^ — log / |IIdgn 



Udq 2 



exp[- / T dt(^ (lq j (t) 2 -e 2 V(q j (t)))+e 2 V R (q 1 (t),q 2 (t)) 



J=l,2 
rT 



J=l,2 



ds 



dt qj (s) ■ W ( qj (s) - qj(t), s - t)qj(t) 



T 



ds / dt qi (s)-W R ( qi (s)-q 2 (t),s-t)q 2 (t) 



(3-1) 



Here J [nd<7 3 -(-)l is the "sum over all paths" and Wr is the photon propagator 

1 



W R (q,t) 



dk\<p(k)\' 



2u{k) 



(t-\k)(k\)e 



-u(k)\t\^ik-r e -ik-q 



(3.2) 



as a 3 x 3 matrix. Here 11 is the unit matrix and \k)(k\ the orthogonal projection 
onto k, k = k/\k\. 

Correspondingly for a single electron 



E 



lim — - log 

T^oo 2T 



ndg(- 



cxp 



dt(\q{t) 2 -e 2 V{q(t))) 



ie 2 

2 C 



ds 



dtq(s)-W {q(s)-q(t),s-t)q(t) 



(3.3) 



Therefore 
£(i?) - 2E 



lim — log E x E 

T^oo 2T 



exp 



— e 



dtV R ( qi (t),q 2 (t)) 



— e 



T 



T 



ds / dtq 1 (s)-W R (q 1 (s)-q 2 (t),s-t)q 2 (t) 



(3.4) 



Here q\{t) and g 2 (t) are two independent copies of the ground state process q(t) 
with a path measure as written in (13.31) . q(t) is stationary and the distribution 
of q(t) at a fixed time is the electronic density computed from the ground state 
ip of H. The average with respect to the ground state process is denoted by E[-] 
and the average over the two independent copies by E x E[-] 
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We note that in the expression (13 .4p the only R dependence sits in Vr and 
Wr, which in a certain sense are small. Thus it is natural to use the cumulant 
expansion. Denoting the exponent by Xt{R), one arrives at 

-^logE x E[e- e2 ^] = -±{?C X {R,T) + \e*C 2 (R,T) + ■■•), 
Ci(i?,T) = — E x E[X T (R)\, 

C 2 (R, T) = E x E[X T (R) 2 } -Ex E[X T (R)} 2 . (3.5) 

As we will see below E x E[X T (R) 2 ]/2T = 0(R~ 6 ), not using the cancellation 
between Vr and Wr terms. Since we are heading for a decay as R~ 7 for large 
R, it should suffice to stop the expansion at the second cumulant C 2 (R, T). This 
will be our main assumption. E x E[Xj<(R)]/2T is exponentially small. The large 
R behavior of the second cumulant will be investigated in detail in the following 
section. 

We note that the expectation under the /c-integral always factorizes with re- 
spect to qx,q 2 - Using their independence we will need to compute only a few 
expectations for the ground state process of a single electron. They are listed now 
for later convenience and proved in Appendix \M We set e — e(k, A), Ej = e(kj, \j) 
for either Xj = 1 or Xj = 2 and 8(t) the step function, 9{t) = —1 for t < 0, 6(t) = 1 
for t > 0. 

(i) E[e ifc -" (<) ] = (^,e ifc -^), 

(ii) E[e ik - q{t) {e{k,X)-q{t))] = -(*fj,e ik - x {H - E){e{k, X) ■ x)tfj) = 0, 
(hi) E[e-' lkl - g{s) e lk2 - q{t) ) = (^j,e- lkl - x e- ls - tKH - E) e lk2 - x ip), 

(iv) E[(ei • q(s))e- [kl - q{s) (l - e^' 2 ^)] 

= 6{t - s)(ip, (ex • x){H - E) e -' lkl - x e- ls - tl{H - E \l - e"^ 2 ^)^), 

(v) E[(ei • g( S ))e- ifcl ^e- ifc2 - g W(e 2 • q(t))} 

= -(V, {e x ■ x)(H - E)e- [kl - x e- ls - tl{H - E) e- [k2 - x (H - E)(e 2 ■ x)tjj) 
+ 5{t - s){e 1 ■ s 2 )(ijj, e- ikl - x e~ ik2 - x ^). 

In (v) the second term arises because locally q(t) is like a standard Brownian 
motion for which dq(s) ® dq(t) = 1.5 (s — t)dsdt. 

The first cumulant can be dealt with immediately. Using (ii) in the above list 
one obtains 

E x E[ qi (s) ■ W R ( qi (s) -q 2 (t),s- t)q 2 (t)] = 0. (3.6) 
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For the potential term it holds 
ExE[V R ( qi (t),q 2 (t))} 

= [ dk\m\Mkr 2 e ik - r ^di-e- ik -^)^d^-^n^) 

JW. 3 

= dx dy((g-5)*if){x)(47r\x + r-y\)- 1 ({g-5)*if)(y). (3.7) 
Jm 3 Jm 3 

Here g(x) is the electron density for ip, which is known to have an exponential 
decay [IS] . 5 is the Dirac delta and * denotes convolution. Since L 3 dx (g(x) — 
S(x)) = and since g decays rapidly, by Newton's theorem it holds that there 
are suitable constants ci,C2 such that 

|E x E[V R ( qi (t),q 2 (t))}\ < Cl e- C2K . (3.8) 

4 The second cumulant 

The second cumulant consists of three terms. They are all proportional to 2T 
by the stationarity of qi(t),q2{t)- We could take the limit T — ► oo first, but the 
symmetric version is more convenient. All expectations are written in terms of H 
and its ground state ip. For notational symplicity we replace H — E by H, hence 
Hip = 0. For inverses as (<pi, H~ n <p2) we make sure that either (<pi,?p) = or 
(02, V') = 0- Note that the ground state is nondegenerate [IE]. But (0i, H~ n (p2) 
could still be infinite. If H is replaced by H at , then H at has a spectral gap and 
therefore an inverse on the orthogonal complement of ip at . 
We set 

C 2 (R, T) = Iyy + 2I VW + I ww (4.1) 

and compute each term separately. The second step is a partial time-integration 
through which one can understand how the R~ 6 decay from the interaction po- 
tential is canceled. In a final step we collect terms according to their number of 
time-integrations and discuss their independence. 
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4.1 Expectations 

a) lyy. Setting 0(kj) = <fj,u(kj) = Uj, we have 

i vv = J dhdt 2 J dM^i^iH^Pd^H^i 2 )^ 1 ^ 

||E[(1 -e-^-" 1 ^ 1 ^! -e- ife2 -" l(i2) )]| 2 



x 



ife2-gi(*2) 



|2 



|E[(1 - e - ifcl - 9l(tl) )]E[(l -e" 

dtxdt 2 J dhdk 2 |£i| 2 |£ 2 | 2 (N%| 2 )-V( fcl+fc2 >- r 

x { |(^,(1 -e- ifcl - x )e- | * 1 -* 2|H (l-e- ifc2 ' a; )^)| 2 

- \ (4j,{l-e- ikl - x )tP)(4j,{l-e-' lk2 - x )iP)\ 2 y (4.2) 

Note that the integrand in (14.21) decays to zero, since 

e-l* 1 -* 2 ^ -> |V)(V>| (4.3) 

as \ti — t 2 \ — > oo. 

b) /yvK- Setting Sj = e(kj, Xj) and noting 

Hk,X))(e(k,X)\ = l-\k)(k\, (4.4) 

A=l,2 

one has 



J yH , = f dhdt 2 dt 3 [ dA; 1 dA; 2 |^ 1 | 2 |^ 2 | 2 (2cu 1 |A; 2 | 2 )- 1 e i ( fel+fc2 )-V 



Ai=l,2 ' 

x E[(ei • gi(ti))e- ifcl - ,?l{tl) (l - e^ 2 ' 91 ^)] 
x E[(ei ■ q 2 (t 2 ))e ikl - q2{t2 \l - e* 2 "^)] 

= J dt!dt 2 dt 3 [ dfc 1 dA; 2 |^ 1 | 2 |^ 2 | 2 (2^ 1 |A; 2 | 2 )- 1 e i(fcl+fc2) - r e- wl|il - i21 

J A 1= l, 2 ^ 

x dfa - t 3 )(i>, (ei • x)iJe- ifcl - x e- |tl " f;i|ff (l - e - ifc2 ' x )^) 

x 6{t 2 - t 3 )(i>, (ei ■ x)He ikl - x e- lt2 ~ t3lH (l - e 1 ^)^). (4.5) 

As proved in Appendix [Bj it holds 

(V>, (ei • x)He~ ikl - x ip) = i{ip, (ei • (p - e^x))^ 1 '^) = 0. (4.6) 
Thus no truncation of the expectation is needed. 
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c) Iww- One has 

W= j dtx&h&h&U f dA; 1 dA; 2 |^ 1 | 2 |^ 2 | 2 (2cu 1 2cu 2 )- 1 e i ( fel+fc2 > r 

x e -^il*i-* 2 | e -^|t3-t4l E J( £i . q 1 (t 1 j)e- ikl - qiitl) (e 2 ■ gi(t 3 ))e~ ife2 " 9l( ' 3) ] 
x E[( ei • q 2 (h))e ik ' q2{t2 \e 2 ■ q 2 (U))e ik ^^] 

dt!dt 2 dt 3 dU J dk l dk 2 \^ l \ 2 \^ 2 \ 2 (2uj l 2u 2 y 1 e< kl+k ^ r 

x e -Wl|fl-t2|g-W 2 |t3-t4| 

x ( - (ei ■ x)He- ikl - x e- ltl -^ H e- ik2 X H{e 2 ■ x)ijj) 

+ 5{t l -h){e l -e 2 )^ ) e-^ +k ^)) 

x ( - (ij>, {e 1 ■ x)He ikl x e- lt2 - ulH e ik2 - x H{e 2 ■ x)ij)) 

+ 6(t 2 - U) ( £l ■ e 2 ) ty, e i ( fel+fe )-^)) . (4.7) 

4.2 Partial time integration 

The next step is a partial time integration for Iyw an d Iww- F° r the integrand 
of Iyw w e use the identity 

6(s - t)(ij, [ex ■ x)He- ikl - x e- ls ~ tlH (l - e^ 2 ^) 

= -%-M>, (e x ■ x)He- ikvX H- 1 e-^ H (l - e - ifc2 ' a » (4.8) 

OS 

and for the integrand of Iww the identity 

- (V>, (ex ■ x)He-' lkl - x e- ls - tlH e- ik2 - x H(e 2 ■ x)ijj) 

= ^ J^, (ei " x)^e- ifcl -^- 2 e-l s -^e- ifc2 (e 2 • 

- 2£(s - t) (ip, (ex ■ x)He- ikl x H- 1 e- ik2 X H(e 2 ■ x)ip) . (4.9) 

We insert these expectations in I VW ,I WW , integrate partially in time, and use 

- u \,-t\ = _ w 2 e -«|«-t| + 2uJ s(s - t). (4.10) 

as at 

The boundary terms vanish. For Iyw one obtains 

Ww= J dtidt 2 dt 3 J] y dA; 1 dA; 2 |^ 1 | 2 |^ 2 | 2 (2 Wl |A; 2 | 2 )- 1 e i(fcl+A;2) ' r 

Ai 

x ( - wfe-^l* 1 "* 3 ' + 2ux5(tx - t 2 )) 

x (i/j, (ex ■ x)He~ ikvx H- 1 e~ ltl ~ t3]H (l - e"^)^) 

x (if>, (ex ■ x)He lkl - x H- 1 e- lt2 - t3lH (l - e 1 ^)^). (4.11) 
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Finally Iww is given by 

Iww = J dtidt 2 d* 3 dt 4 J &h&k 2 \ip 1 \ 2 \$ 2 \\2u l 2uj 2 y 1 e^ +k * ) - r 

X { ( - w^^ltx-tal + 2ui5 ( ti _ h ^ ( _ ^2 e -. 2 |t 3 -t 4 | + _ t4 )) 

x (ei • x)He- ihl - x H- 2 e~ ltl - t3lH e-' lk2 X H{e 2 ■ x)ip) 
x (ei • x)He ikl - x H- 2 e-^- ulH e ik2 - x H{e 2 ■ x)ip) 
+ 0^ - i 3 )wi0(t 3 - t 4 )^2e"' JJl|tl ~* 2| e~" 2|i3 - t41 

x (ei • x)He-' lkl X H- 2 e- ltl - t3lH e- ik2 - x H{e 2 ■ x)ip) 

x <J(t 2 -U)B + (k u Ai;fc 2 ,A 2 ) 

+ (£i • x)He ikl X H- 2 e- lt2 - tilH e ik2 - x H(e 2 ■ x)^) 

x 5(h - t 3 )B4h, Ai; fc 2 , A 2 )) 

+ e -^|tl-i2l e -^2|t3-t4| 5 ( ti _ _ ^) 

x Ai; k 2 , X 2 )B+{k 1 , A i; k 2 , A 2 )}, (4.12) 

where 
B-(ki, Ai; fc 2 , A 2 ) 

= ( £l •e 2 )(^,e- i(fcl+fc2) » -2(V',(ei-x)i7e- ifcl - x i/- 1 e- ifc2 ' :c iy(e 2 -x)^), (4.13) 
5+(A;i, Ai; fc 2) A 2 ) 

= (ei • e 2 )(i;,e i{kl+k2) - x ip) - 2(^, ( £l • x)He ikvx H~ 1 e ik2X H(e 2 ■ x)if>) (4.14) 

We use (14.11) and collect the terms of C 2 (R, T) according to the number of their 
time-integrations, divide by 1/2T, and take the limit as T — > oo. To prepare for 
the limit R — > oo, we rescale the momentum integration as kj ~» kj/R,j = 1,2, 
and introduce the unit vector h = r/R = (0, 0, 1). Note that e(k/R, A) = e(k, A). 
The two-time, three-time, and four-time integrations are treated separately. 
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4.3 Two-time integrations 

The sum of all terms involving two-time integrations is denoted by I 2 (R). One 
has 



h{R) = I dhdhmh/R^ih/R)^^ 

dt{i?" 2 (|A; 1 | 2 |A; 2 | 2 )" 1 f|^,(l-e- ifcl -^)e-'^(l-e- ifc2 - a; / ii 



- | (rjj, (1 - e" ifc ^/» (if,, (1 - e- ife ^)^> T) 
+ 2 iT 4 N~ 2 (^, (ei ■ x)He~ ikl - x/R H~ 1 e~ WH (l - e^ 2 '*^) 

Ai 

+ ^ /T 6 (^, (ei • x)He- ikl - x/R H- 2 e- m e- [k2 x/R H(e 2 ■ x)ip) 

Ai,A2 

x (ei • x)He ikl - x/R H- 2 e- m e ik2 - x/R H{e 2 ■ x)ij)) 
+ ^ fr 3 (2^ 2 (wi + W2))" 1 

Ai,A2 

x B4h/R, An fc 2 /i2, X 2 )B+(k 1 /R, An fc 2 /i2, A 2 ) 
=h,i(R)+h,2(R)- (4-15) 

We consider the sum, I 2j i(R), of the first three terms and expand in 1/-R, 
which yields expectations of the form (^>, (a • x)e~'''^(6 ■ x)^) with a, 6 G M 3 . By 
rotation invariance 

(V>, (a ■ x)e- |t|H (6 • x)^) = (a • x ■ e" 1 * 1 ^). (4.16) 

Using 

£ ( £ 1 ■ = ! + (*!• ^) 2 , $>l ■ ^ = 1 " & • ^ 417 ) 

Ai,A2 Ai 

one arrives at the lowest order 

R~\\^,x- e-l*l^)) 2 ((fc 1 • k 2 ) 2 + 2 - 2ft • k 2 ) 2 + l + k 2 f). (4.18) 

The nonsmooth contributions, containing (ki ■ k 2 ) 2 , are canceled exactly, while 
the smooth part, at this order, is given by 

hx{R) 

= 3jdtJ dfc 1 dfc 2 |^(fc 1 /i2)| 2 |^(fc 2 /^)|V (fcl+fc2) ".R- 6 (|{^,a: • e" 1 * 1 ^^)) 2 , 

(4.19) 
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which inherits the rapid decay in R from ip. At the next order one picks up the 
quadratic contributions kf, ki ■ k 2 , k\ with coefficients integrable in t. By power 
counting one arrives at I 2j i(R) — R~ 8 as R — > oo. 

For the second summand, I 2y2 (R), we use that to leading order in 1/R, 

B^h/R^X^h/R, A 2 ) = ^I^H^Ke! • k 2 ){e 2 ■ h)^^ 2 ^)) + 0(iT 4 ), 

(4.20) 

as in shown in Appendix IB1 Using that |y3(0)| 2 = {2ir)~ 3 one arrives at 

= ^ 7 (|(^,xV)) 2 (27r)- 6 i J dk 1 dk 2 {u 1 u 2 {u l + u^)' 1 ^ 2 ^ 2 
x e i(fcl+fc2)A (l - (ki ■ kff + 0(R~ 9 ) 

= — (±-) 2 e-*alR- 7 + 0(R-»). (4.21) 

IX 111 

The evaluation of the numerical coefficient is discussed in AppendixO Altogether 

hmi? 7 J 2 (i?) = ^(^)V 4 ^. (4.22) 

4.4 Three-time integrations 

The sum of all terms involving three-time integrations is denoted by ^(-R). One 
obtains 

2 \,n!b„l EA|2_i(fci+fc 2 )-n 



h{R)= J dhdhj dhdhlfiih/R^lvikz/R^e 
[-R'^lk^e-^-^ 



x 
x 



(y^ty, (ei • x)He- [kl - x/R H- 1 e- ltllH (l - e"^ 2 ^)^) 



x (ip, {e x ■ x)He ikl x/R H- 1 e- lt2lH (l - e ik2 ' a 
- R-Wu^e-^-^* + \u 2 e-^-^ R ) 
x ( ^ (ei • x)^-^-^^-^- 1 * 11 ^-^^//^ ■ x)V) 

Ai,A2 

x (ei • x)He ikvx/R H- 2 e-^ H e ik2 - x/R H(e 2 ■ x)t/j) 
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+ J R~ 6 ^(t 1 )e" wl| * l|/ ^(t 2 )e- W2|i2|/R 

x ((^,(e 1 -x)He- ikvx/R H- 2 e-\ tl - t2 \ H e- ik2X/R H(e 2 -x)i>) 

Ai,A2 

x B + (A;i, Ai;fc 2 ,A 2 ) 

+ (V, (ei • x)He ikl - x/R H- 2 e- ltl - t2lH e ik2 x/R H{e 2 ■ x)ip) 
x Ax; fc 2 ,A 2 )))} 

=I 3A (R)+I 3t2 (R). (4.23) 

We expand the integrand in which yields 

d*id* 2 y dfc 1 dfc 2 |^(fc 1 / J R)| 2 |^(A; 2 / J R)| 2 e i ( fcl+fc2 >" 

x { - rVe" ui|fi " bl/fi (i - (ifei ■ 4) 2 ) 



x x ■ e- lhlH xip)l(ip, x ■ e- lt2lH xip) 
- R- 7 {\u l e^\ t ^/ R + \u 2 e-^-^' R ){\ + (h ■ k 2 ) 2 ) 
X x ■ e-^ H xt/j)l(i/j, x ■ e~ lt2lH xip) 

+ R- 6 \e{t 1 )e-^ t ^ R e{t 2 )e~ uJ2 \ t2 \i R 

x ( K^^-e" 1 * 1 "* 2 '^) 

Ai,A2 

x (e 1 ■ e 2 ) (B+ih/R, A l5 k 2 /R, A 2 ) + A l5 k 2 /R, A 2 )) 

= -R- 7 J dhdhj dA; 1 dA; 2 |^(A; 1 /i?)| 2 |^(A; 2 / J R)|V (fcl+fc2) -Ve- aJl|tl -* 2|/ii 

+ i?- 6 y dhdt 2 J dhdhlipih/R^lrth/R)^^ 

Ai,A2 

x e i(fcl+fc2) -"^(ti)e- Wl| ' l|/ ^(t 2 )e-" 2| ' 2|/i? |(^, x ■ e~ ltl ~ t2lH xtP) 

x ( £l ■ e 2 ) (B + (h/R, Ai; fc 2 /i?, A 2 ) + B_{h/R, A i; fc 2 /#, A 2 ))) . (4.24) 

For the first summand the nonsmooth contributions are canceled exactly, as in 
Section l4~3l while the smooth contribtions decay rapidly since ip does so. Its next 
order picks up an extra factor R~ 2 . Thus I 3t i(R) ~ R~ 9 for large R. For the 
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second summand we perform explicitly the time integration with the result 

h t2 (R) =R~ 6 J dhdk^kjR^ih/R)^^^ 

x R(uj x + ^ 2 )- 1 |(^, x ■ (H(H + irVr 1 ^ + RT 1 ^)- 1 )^) 
x ( £l ■ £2 ) ( B +( k i/ R > A i! A 2 ) + B_(h/R, Ai; fc 2 /fl, A 2 )) 



(4.25) 

Since 5± = -R~ 2 |/ci||fc 2 |(ei-fc 2 )(e 2 -A;i)(|('^,x 2 '^))+C(i?~ 4 ), as shown in Appendix 
iBl we conclude that 

h >2 (R) = - ir 7 (§(^,arV)) H(i/>,x ■ H- 1 ^)) 

x 2(2ir)- 6 J dA;idA; 2 e i(fcl+fc2 >"(a;ia; 2 (a;i + a; 2 ))- 1 
x |fci| 2 |fc 2 | 2 (fci • fc 2 )(l - (fci • £ 2 ) 2 ) + 0(iT 9 ) 
=— (t^) V 4 a E a M iT 7 + 0(fT 9 ). (4.26) 

The evaluation of the numerical coefficient is discussed in AppendixO Altogether 
one has 

lim R 7 I 3 (R) = — (^) 2 e- 4 Q E a M . (4.27) 

R— too 7T 271 

4.5 Four- time integrations 

There is only a single term with four-time integrations, namely 



h{R) = J dtidt 2 dt 3 / d^id^2 j ! 2 1 ^2 1 2 e i(fel " 



■k 2 )-r 



Ai,A2 

x - Wl w 2 e" Wl| * 1 - ta+ * s| e -W3| * s| 

x (ei ■ x)He- ikl - x H- 2 e- lhlH e- ik2 X H(e 2 ■ x)rj)) 

x {e x ■ x)He ikl - x H- 2 e- lt2lH e ik2 X H(e 2 ■ x)i>). (4.28) 



We rescale kj ~> kj/R as before and in addition t$ ~> t$R. Then the exponentially 
decaying terms are 

e -Wl|t3+(tl-t2)/fl| e -W2|t3| ^ g-Wl|*3lg-W2|*3| (4.29) 
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for large R. Expanding in 1/R yields 

dt / dk 1 dk 2 \ip(k 1 /R)\ 2 \ip(k 2 /R)\ 2 

x e i(fcl+fc2) -V^ 2 e"* K+W2) (l + (h ■ hf) + 0{FT 9 ) 
=R' 7 (2ny 6 (l(ij,x ■ H~ 1 xi))) 2 2 J dk^k 2 
x el():i+fc2 ). A J^ L(1 + ^ . ^ )2) + 0(jR _ 9) 

CUi + UJ 2 

= R ~ 7 ^ 2e ~ 42 i a - +0{R ' 9) - (430) 

The prefactor of R~ 7 agrees with the one computed already by Casimir and 
Polder. 

4.6 Sum of all terms 

We add the limits listed in fH~22l) . (Q7I> . and fOOll . which yields 



K ( e ) = ^ fl7 (^^ 4 i C2 < fl - r ») = (i) 2 ( 



23 2 26 64 2 

— — «E H «E«M H «M 

167T 7T 7T 



(4.31] 



4.7 Comparison with previous results 



Magnetic contributions to the —R~ 7 decay were first considered by Feinberg and 
Sucher [SJ, [9] and by Boyer [10] . They find that a| at and at have the same 
coefficient, namely 23/167T, while the one of «E,atttM,at is 7/8n. This raises the 
issue on the origin for the discrepancy. 

The q| at term can be most easily obtained through the dipole approximated 
hamiltonian 

#dip =\{pi - eA(0)) 2 - e 2 V{x x ) + ^(p 2 - eA(r)) 2 - e 2 V(x 2 ) + H t 

+ e 2 V R (x 1 ,x 2 ). (4.32) 

One shifts p\ by eA(0) and p 2 by eA(r) through the unitary U = exp[— i(xi • 
A(0) +x 2 ■ A(r))]. Then 

UH^U' 1 = l -p\ - e 2 V( Xl ) + l -p\ - e 2 V(x 2 ) + H t 



+ e 2 



exi ■ E±(0) - ex 2 ■ E±(r) + e 2 V R (xi, x 2 ) 

j dk\m\ 2 (^A + \x 2 2 +^ k - r {x l -x 2 - \k\- 2 ( Xl ■ k)(x 2 ■ *;))), 

(4.33) 
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where Ej_(r) is the quantized transverse electric field. Note that the long range 
part of Vr is cancelled. From the 4-th order perturbation in — ex\ ■ E±(0) — ex 2 ■ 
E±(r) one obtains the Casimir- Polder result. 

To include magnetic effects one proceeds to the next order of the multipole 
expansion and defines 

#mui =\{pi - eA{0) - ex! ■ V r A(0)) 2 - eV(zi) 

1 2 

+ - (p2 - eA(r) - ex 2 ■ V r A(r)) - e 2 V(x 2 ) 

+ H t + e 2 V R {x 1 ,x 2 ), (4.34) 

which, as before, is unitarily transformed to 

UHrnaU- 1 = l -( Pl - ex, ■ V r A(0)) 2 - e 2 V(x x ) 

1 2 

+ g m - ex 2 • V r A(r)) - eV(x 2 ) 

+ tff - ex! ■ E ± (0) - ex 2 ■ E ± (r) + e 2 V%*(x u x 2 ), (4.35) 

where V R T is the interaction potential from (14.331) . Note that V x A (x ■ S7 r A(r)) = 
B(r) = V r A A(r). Thus the hamiltonian (I4.35P clearly displays the fluctuating 
electric and magnetic fields. Expanding in the terms proportional to V r A(0), 
V r A(r), E±(0), E±(r) to 4-th order leads to an energy with a large R asymptotics 
in agreement with ( 14.31H . In spirit the authors of [8], El [10] make a further gauge 
transformation through the unitary exp[— i(xi ■ V r (xi ■ A(0)) + x 2 ■ V r (x 2 • A(r)))]. 
Thereby X\ ■ V r A(0) is transformed to A B(0) and correspondingly for x 2 ■ 
V r A(0). In addition there are terms coming from the shifting of E±. In the 4-th 
order perturbation only the former terms are taken into account. More precisely 
the term (0 a t ® 0, {xi A B(0)) 2 PH~ 1 P(x 2 A B(r)) 2 4> at <g> Q), which yields indeed 
the integrand ( 1C.1I) and thus 23/167T for the prefactor. Here P = 1 — \ip)(ip\ and 
Q is the Fock vacuum. For the cross-term they use (ip at <g>Q, x± ■ E_ s _(0)PH~ 1 Pxi ■ 
E ± (0)PH- 1 P(x 2 A B(r)) 2 %l). At ® Q), which yields the integral d2l) with the term 
(1 — (k\ ■ k 2 ) 2 ) omitted and thus 7/8n for the prefactor of aE,at«M,at- 

5 The dimensionless strength 

We restore the physical units in H of (12.11) . Then 

H = ^(P- ~c A &)) 2 ~ e 2 V(x) + ^ (5.1) 

with p = —iftV x , 

H f =J2[ dkhc\k\a{k,\)*a{k,X), (5.2) 
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and 



A(x)=J2 [ dk^/^k\^(k)e(k,\)(e ik - x a(k,\)+e- ik - x a(k,\y). (5.3) 

A=1,2- /r3 

H is transformed to atomic units through the canonical transformation U defined 

as 

U*a(k, X)U = {a- 2 X c ) 3 / 2 a{a- 2 X c k, A), 

U*xU = a~ 1 r B x, U*pU = ar^p. (5.4) 

Then 

U*HU = a 2 mc 2 (^(- iV x - V^a 3/2 A(ax)) 2 - V(x) + H^j 

= a 2 mc 2 H, (5.5) 

where 

Hf=V] I dk\k\a(k,X)*a(k,X) (5.6) 

A=l,2 J* 3 

and 

V(x)=4n [ dk^X^kflkl^e-' 1 ^, (5.7) 

A(x) =Y [ dk ^ a ^l_ k K (k, X) (e ik ' x a(k, X) + e~ ik - x a(k, A)*) . (5.8) 
A T7/ 2 Ji83 \/2\k\ 



Energy is now in units of a 2 mc 2 and distances are in units of the Bohr radius r-g. 
From (15. 6p . (I5.8P the photon propagator in atomic units is obtained as 

W R (q,t) = 4* f dfc|^(a 2 A c - 1 fc)| 2 ^(]l-|fc)(fc|)e-"" 1 l fc ll*W fcr e- i ^. (5.9) 
Jr 3 ^l fc l 

Comparing with (13.21) this amounts to replacing e 2 by Att and t by t/a. One 
merely has to follow through this change in the computation of Section HI The 
final result is 

E(R) -2E = -K{a, X c A)amc 2 (R/r B y 7 (5.10) 
for large R with the strength k obtained as 

k(cx, A c A) = - ( ^<Se + 26<5 E a M + 645^) . (5.11) 

7T \ lb / 

Here, the dimensionless electric and magnetic dipole moments are 

ti E = 2(l(i),x- (H - E^xip)) (5.12) 
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and 



a M = -a 2 -(|(^,x 2 ^)) 



(5.13) 



with ijj the ground state of H, Hip = Eip. The ultraviolet cut-off is implemented 
by replacing 0(k) by the form factor <pi(k/A): 1 decreases rapidly at \k\ = 1 
from (2tt)~ 3 / 2 to 0. 

Following Bethe [19] a physically reasonable choice is A C A = 1. One would 
like to remove the ultraviolet cut-off through A — > oo. But this limit is not well 
understood. In any case, the bare mass m would have to substituted by the 
renormalized mass. For fixed A, A C A = 1, in the limit a — > 0, H decouples 
from the radiation field and n(a, A C A) tends to the strength obtained by Casimir 
and Polder. Then qe ~ 9/4 and <5m — —& 2 with quality for the strict Coulomb 
potential. For a systematic expansion in a one would need the correction of order 
a to (5e and to the bare mass m. 

6 Conclusions and Outlook 

The reader may wonder how much is still missing to a complete proof. For the 
full model, a central difficulty is that the ground state process q(t) is not so well 
under control. Since H has no spectral gap, the two-point correlation has a slow 
decay, presumably as K[q(t)-q(0)} ~ |t|~ 4 for large t. This means that higher order 
cumulants are difficult to control. In fact, we cannot even prove that qe < oo. 

From the statistical mechanics point of view an interesting case would be to 
replace the ground state process q(t) by g a t(£), namely the one governed by H at . 
This process is exponentially mixing which should help in controlling the error 
term. On the other hand, the double stochastic integral in the action causes extra 
difficulties. Unfortunately there is no obvious hamiltonian corresponding to this 
approximation. 

A further variant would be the dipole approximation of equation ( 14.32ft . In 
the path integral (13. 4p this corresponds to replacing Wr(x, t) by W(0,t). The 
effective action is quadratic in <jj(t) and the partial time integration can be easily 
implemented. It results in a diagonal term which cancels the slowly decaying part 
of Vr and the remainder action is given by J dsdtq(s) ■ W R (0, s - t)q(t). Since 
Wr{Q,£) ~ t~ A for large t, the difficulties mentioned above remain. To have an 
exactly solvable Gaussian model, on top one would have to use the quadratic 
approximation for V and Vr. 

For ground state properties a powerful method is the Feshbach projection 
together with a successive integration over high /c-modes of the radiation field 
[201 |2T] . It would be interesting to understand whether this technique could be 
used for a rigorous control on the van der Waals potential. 

Acknowledgement: We are grateful to P. Milonni for very useful hints on the 
literature. 
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A Proof of (i) to (v) 

The proof of (i) to (v) is based on the identity 

m 

®[HfMtj))] =^Ji(^)e H4l - 42|(H - S) ---/ m -i(x)e- | ''™- 1 - < - |(H - i?) / m (x)^) 

(A.l) 

for the time order t\ < t 2 < • ■ ■ < t m . q(t)dt is the Ito stochastic integral 
as defined through the forward discretization. For this purpose we introduce 
the lattice spacing 5 and let [t]g be the integer part of t modulo S. Then 
q(t) = lims-,.0 5~ 1 (q([t]s + 5) — q([t]s))- We only establish (ii) and (v). The 
other items are proved by the same procedure, 
ad (ii): Using stationarity, 

E[e ik - q{t) e ■ q{t)\ = lim ^E[e ihq ^e ■ (q(5) - q{0))\ 

= lim h^, e ^u-s(H-E) _!)( e . x u\ 
= -(^e ik - x (H-E)(e-x)^), (A.2) 
which vanishes as proved in Proposition IB. 41 

ad (v): For t ^ s we proceed as in ad (ii), which yields the first summand of (v). 
For the "diagonal part" one has, for small 5, 

E[( 6l • (q(5) - q(0)))e- ik ^e- ik ^°\e 2 ■ (q(S) - q(0)))] 
= {^^- ik ^ e - ik ^{e-^ H - E )( £l . z )( e2 . x) - (ei • x)e~ S{H - E \e 2 ■ x) 
- (e 1 ■ x)e- 5 ^ H - E \e 2 ■ x) + ( £l • x)(e 2 ■ x)e- 5 ^ H ~ E ^) 
= -5(^,e- [kl - x e- ik2 - x [[5~\l - e~ 5{H ~ E) ), {e x ■ x)}, (e 2 • x)]^j) 
~5(e 1 -e 2 ){iP,e- ikl - x e- ik2 - x iP). (A.3) 

B Ground state expectations 

In this appendix, we mainly work in the Schrodinger representation L 2 (K^) ®$ = 
L 2 (Rl) ® L 2 (Q), see [TH [TB]. In this representation it holds 

(a) A(x) is a real- valued multiplication operator. 

(b) Let $ = r(e 17r//2 ). Then ^e - *^?? -1 is positivity improving. Hence ip = #0 is 
strictly positive. 

Here for a unitary operator U, T(U) is defined by 

r(c/)a(/o* • • • a(f n yn = aiuhy ■ ■ ■ a(uf n yn 
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for all n G N and /i, . . . , f n G C 2 (8)L 2 (M 3 ), where Q is the Fock vacuum. A linear 
operator A "improves the positivity" if, for all < / G (L 2 (Rl) ® L 2 (Q))\{0}, 
e -tAy > q for gj_j > q L e t j be the natural conjugation in L 2 (Rl) ® L 2 (Q), 
namely, = \l/ for all ^ G -^ 2 (^.) <8> L 2 (Q). A linear operator A is called to be 
J-real if A commutes with J, i.e., J A = AJ. Then, from (b), it follows that 

(c) H = r dH r d~ l is a J-real operator, i.e., JH = HJ. 



Lemma B.l (Vanishing ground state expectation I). Let F(x) a measurable func- 
tion on IR 3 such that ip G dom(F). Then one has 

{tp, F{x)A{x)^) = 0. (B.l) 

Remark B.2 From the proof one infers the stronger property 

F(x) A(x) n e-^ H ■ ■ ■ e- s ^ H A(x) t2n+1 ?P) = 

for all n G {0} U N, s x , ■ ■ ■ , s 2n > and F G L°°(M 3 ). (Note that we have to use 
(c)0 

Proof. Note that 

M = r(e i7r )tf J. (B.2) 

Since A{x) is a real- valued multiplication operator, one has A(x)J = JA(x). 
Hence using (IB. 21) . one sees that E(x) = r dA(x)'d~ l is purely J-imaginary, that 
is, JE(x) = —E(x)J. Remark that 

ftM(a:)V) = $,E(x)$). 

Since ip is real- valued by (b), the right hand side is purely imaginary. On the 
other hand, by the self-adjointness of A(x), the left hand side is real. Thus the 
only possibility is (ip, A(x)tp) = 0. Similarly one can show that 

(iJ;,3tF(x)A(x)il>) = = (tp, QF(x)A(x)ip) 

which implies the assertion. □ 

Lemma B.3 (Vanishing ground state expectation II). One has the following: 

i) For all a G M 3 , (ip, (a ■ x) 2n+1 tp) = 0. 

ii) If m + n is odd, then (ip, (a ■ x)H(b ■ x) m H -1 (c • x) n H(d ■ x)tp) = for all 
a,b, c,d G IR 3 . 
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Proof. For the proof, we return to the Fock representation L (R^) Let J 2 be 
the involution defined by 

J 2 ^=V * n (-2c;fci,Ai,...,A^,A n ) (B.3) 



for all * = J2n> ^n(x; fa, Ai, . . . , A; n , A n ) G L 2 ^,) ® Then as proved in [22], 
we can check that 

J 2 x = -xJ 2 , (B.4) 
J 2 H = HJ 2 , (B.5) 

namely, x is purely J 2 -imaginary and H is J2-real. Note that J 2 ip = e lB ip for some 
9 G [0, 27r) by ( IB. 51) and the uniqueness of the ground state. Hence 

(j>, (a ■ x) 2 " + V) = - (a • x) 2n+1 J 2 ^) 

= -(^,(a-x) 2 " + V), (B.6) 

which implies (a • x) 2n+1 ^) = 0. Similarly if m + n is odd, then 

(V, (a • x)H(b ■ x) m H^(c ■ x) n H(d ■ x)i/j) 

= -{J 2 iJ>, (d ■ x)H(c ■ x) m H~\b ■ x) n H(a ■ x)J 2 tfj) 

= -ty, (d ■ x)H{c ■ x) n H-\b ■ x) m H(a ■ x)i/>). (B.7) 

On the other hand, in the Schrodinger representation, we can see that 

(d ■ x)H{c ■ x) n H-\b ■ x) m H{a ■ x)$) 
= (a ■ x)H{b ■ x) m H^(c ■ x) n H(d ■ x)$) (B.8) 

because every operator appearing in the expectation is J-real. Therefore 

{i/j, (d ■ x)H(c ■ x) n H-\b ■ x) m H(a ■ x)ip) 
= (tfj, {a ■ x)H(b ■ x) m H-\c ■ x) n H{d ■ x)ip). (B.9) 

Combining this with flB.Tj) . we conclude ii). □ 

Proposition B.4 One has the following 

i) (4j,e ik - x H{e{k,X) -x)iff) =0, 

ii) B±{fa/R,\i,k 2 /R,\ 2 ) = RT^faWH^ ■ k 2 )(e 2 ■ fci)(|(V,^» + O(R^) 
as R — > 00. 
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Proof, ad i) By (1B.13P below and by (ip,e lk ' x A(x)i/j) = according to Lemma 
IB.1[ one has 

(iff, e ik ' x H(e(k, A) • x)ifj) = - i(ip, e ik ' x (e(k, A) ■ (p - eA(x)))0) 

= -i^ ) e*-(e(fc ) A).p)^> 

= - ~J dxdfi(A)e ik ' x e(k, A) • (-iV^ 2 )(x, A) 

= ~\j dxdfi(A)e lk - x (e(k, A) • £;)^ 2 (x, A) 
=0, 

since • s(k, A) = 0. 

ad ii) By Lemma [B. 3 1 the order vanishes and it suffices to check the order 
R° and R~ 2 . 

ii-a) order R°. One has 
lim B.ih/R, Ai; fc 2 /#, A 2 ) = (ei • e 2 ) - 2(0, (ei • x)H{e 2 ■ x)ijj). (B.10) 

Furthermore, using the identities 

[a ■ p, b ■ x] = —ia ■ b, (B.ll) 
[H,a-x] = -ia-p, (B.12) 
H(a ■ x)ip = — ia ■ pip (B.13) 

for all a, b G C 3 , where p = p — eA(x), one obtains 

2(0, (ei • x)H(e 2 ■ x)ip) = i(ip, ((e x ■ p)(e 2 ■ x) - (e x ■ x)(e 2 ■ p))ip) = (e 1 ■ e 2 ). 

(B.14) 



Remark that we have used ( IB. II) to conclude that (0, (ei • y4(x))(e 2 • x)0>) = 
= (ip,(e 1 -x)(e 2 - A(x))iP). 

ii-b) order i? -2 . We expand as 

R 2 B4h/R, Ai; fc 2 /#, A 2 ) = Cl A; 2 + c 2 £; 2 + d(k u A x ; fc 2 , A 2 ) + 0(iT 2 ). (B.15) 
For Ci/c 2 , the expansion coefficient is 

(e 1 ■ e 2 )(i/j, (h ■ x) 2 ifj) - 2(0, (e 1 ■ x)H(k 1 ■ x) 2 (e 2 ■ x)iP). (B.16) 

Using (IB . ID and that ip and H can be chosen as real, for the second term one 
obtains 

(ip, (ei ■ x)H{ki ■ x) 2 (e 2 ■ x)ip) + (ip, (e 2 • x)(k x ■ x) 2 H(e 1 ■ x)ip) 
= i(0, (ex ■ p)(h ■ x) 2 (e 2 ■ x)ip) - i(ip, (e 2 ■ x)(k x ■ x) 2 (e 1 ■ p)ip) 
= i(ip, (k 1 -x) 2 [(e 1 -p),(e 2 -x)]ip) 

= (e 1 -e 2 )(iP,(k 1 -x) 2 iP). (B.17) 
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Thus C\ = = c 2 . 

The expansion coefficient of the mixed term dfa, Ai; k 2 , X 2 ) is given by 

-(ei ■ e 2 ){ip, (h ■ x)(k 2 ■ x)ip) + 2(ip, fa • ■ x)H~ l {k 2 ■ x)H(e 2 ■ x)ij>). 

(B.18) 



Using fIB.ip and that H and ip can be chosen to be real, we have 



(V>, (ei ■ z)if(fci • x)H-\k 2 ■ x)H(e 2 ■ x)ijj) + (1 <-> 2) 
1 



(£1 • x) (h ■ x) (k 2 ■ x) (e 2 ■ p)i/j) 

- (if;, (ei ■ x)(/ci ■ p)H~ 1 (k 2 ■ x)fa ■ p)ij)) 
+ i(tp, fa ■ p)(Ai ■ x)(fc 2 • ar)(e 2 ■ x)ip) 

+ (V, (£1 • p){k\ ■ x)H~\k 2 ■ p)(e 2 ■ x)i)) + (1^2) 

= ^{^^(h- x)((e 1 -p)fa-x) - fa ■ x) fa ■ p)){k 2 ■ x)ip) 

- (4>, fa ■ x)(ki ■ p)H~ 1 (k 2 ■ x)fa ■ p)^) 

+ (iP, fa ■ p)(k 1 ■ x)H~ 1 fa ■ x)(k 2 ■ p)^> + (1^2)} 

= fa ■ e 2 )(ifj, (hi ■ x)(k 2 ■ x)ifj) + fa ■ k 2 )fa ■ ki)(l(i/j,x 2 ip)) 

= i(^,x 2 V)(( £l • e 2 )(h ■ k 2 ) + fa ■ k 2 )fa ■ h)). (B.19) 

This proves the assertion. □ 

C Numerical coefficients 

In this appendix, we will explain how to compute the following integrals appearing 
in the main text: 

ft = f dhdk 2 l } hllh } e«fci+fa)-*(i + (ki ■ hf), (C.l) 
J \ki\ + \k 2 \ 

S 2 = ( dhdk 2 .j hl } h j . e^+^Ch ■ fc 2 )(l - (ki ■ h) 2 ), (C.2) 
J \ki\ + \k 2 \ 

s *= [ dk ^ jrrM 1 e i{fcl+fc2) "(i - (jfei • hff. (c.3) 

J \ki\ + \k 2 \ 

These integrals are of the form 

Sj = [ dfrdfc j^j . e^+^F^h ■ k 2 ), (C.4) 

J \Ki \ + \K 2 \ 



which we rewrite as 



(C.5) 
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Let us switch to polar coordinates (r, <p, $) by 

k = (Y cos tp, Y sin tp,X), 

X = cos$, F = sin^. (C.6) 
Clearly X 2 + Y 2 = 1. Then we have 

k\-k 2 = cos(^! - ^2)^2 + XlX 2 , (C.7) 

and hence 

^ = / dt / drtrfe-^ 1 / dX a e iriXl 
Jo Jo J-i 

/oo /» 1 

dr 2 r 3 2 e" tr2 J dX 2 ^ X2 & j (X 1 ,X 2 ), (C.8) 



where 



& j {X 1 ,X 2 )= / d^ / dip 2 F j (cos(<p 1 -ip 2 )Y 1 Y 2 + X 1 X 2 ) 
Jo Jo 

r2ir 

= 2vr / di P F j (cosipY 1 Y 2 + X 1 X 2 ). (C.9) 



Set 



»i 

(A) = 2n I drr 3 e" ir / dXe irX A(X). (CIO) 



1 



After performing (^-integration in (10.91) . Sj can be expressed as 

S 1 = jHdf{§<l)<l) - (1)(X 2 ) + |(X 2 )(X 2 )}, (C.ll) 

S 2 =J™dt{ - \{X){X) + 3(X)(X 3 ) - |(X 3 )(X 3 )}, (C.12) 

S 3 = jHdf{§(l)<l) + |(1)(X 2 ) + |(X 2 )(X 2 ) + f <1><X 4 > - f (X 2 )(X 2 ) 

+ f(X 4 )(X 4 >}. (C.13) 



Using Mathematica, one obtains 

—4 + 12/ 

M = W (C - 14) 



(X 3 )^4i{ ^±^^ -3arccot(t)}, (C.17) 

, 4(3 + 27t 2 + 32t 4 + 12t 6 ) A . 

(X 4 ) = -i , — '- - 48t arccot(t). (C.18) 

(l + t i y 
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Inserting these formulas to (10. 1 lH - flO. 13[) and using Mathematica again, one ar- 
rives at 



Si = 92vr 3 , S 2 = 208vr 3 , S 3 = 256vr 3 . (C.19) 
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